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1. 


Recent  Advances  in  The  Edge-Function  Method 


Summary 


Major  advances  in  the  Edge-Function  Method  in  1979-80  were 
made  in  four  main  areas: 

1.  Anisotropic  Thin  Plates  (J.J.  Grannell  and  P.M.  Quinlan). 

2.  Free  Vibrations  of  Thin  Plates -(Richard  Sheehy). 

3.  Analysis  of  Elastic  Plates  and  Bridge  Slabs  (A. A.  Ahmud). 

4.  Elliptical  Cracks  in  a  Prismoidal  Body  (P.M.  Quinlan, 

J.J.  Grannell,  A.N.  Atluri,  J.E.  Fitzgerald). 


The  principal  investigator  was  actively  involved  in  an  advisory 
role  with  both  R.  Sheehy  and  A. A.  Ahmed  on  their  Ph.D.  theses.  The 
present  report  gives  an  abstract  and  conclusions  in  each  of  the  above 
works.  It  also  gives  some  examples  from  (2)  and  (4)  to  illustrate  the 
range  of  E  F.M.  in  Free  Vibrations  and  in  buried  Elliptical  Cracks  in 
Fracture  Mechanics. 


These  examples  are  the  touchstone  by  which  the  present  state  of 
development  and  utility  of  E.F.M.  can  be  judged.  In  all  cases  root  mean 
square  boundary  residuals  are  calculated  as  a  routine,  thus  providing  a 
practical  measure  for  judging  the  acceptability  of  the  solution  proferred. 
Tf  the  residuals  are  within  the  limits  within  which  an  engineer  can  specify 
the  boundary  conditions  of  the  problem,  then  the  corresponding  Mathematical 
Model  which  E.F.M.  offers  is  as  "exact"  a  solution  to  the  physical  model, 
or  problem,  as  can  be  obtained. 


1 
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Grannell  J.J.  and  Quinlan  P.M. 

"The  Edge-Function  Method  for  Thin  Anisotropic  Plate  Bonding" 
Vol.  80A  No.  1  Proc.  Roy.  Irish  Academy  1980 

Abstract 


The  edge-function  method  is  developed  for  the  bending  of  thin 
anisotropic  polygonal  plates.  The  solution  is  represented  as  a 
superposition  of  asymptotic  solutions  of  the  plate  bending  equation  in 
the  neighbourhood  of  characteristic  segments  of  the  boundary.  In  the 
case  of  a  polygonal  plate  these  segments  are  edges  and  vertices. 

Interior  polynomial  solutions  are  included  and  transverse  loadings  are 
modelled  using  a  Green’s  function  approach.  The  boundary  conditions 
are  imposed  using  accelerated  discrete  least-squares  Fourier  interpol¬ 
ation  over  each  of  the  edges.  The  accuracy  of  the  model  is  assessed  by 
computing  the  boundary  resid:  Is.  A  symmetry  principle  is  developed  for 
solving  symmetric  boundary  value  problems.  The  solution  of  three  plate 
problems  is  presented  together  with  details  of  accuracy  and  convergence. 


Three  Plate  Problems  studied 

1.  Orthotropic  Cantilever  Plate  under  uniform  normal  load 

2.  Orthotropic  Skew  Plate  under  uniform  normal  load.  Plate  fully  fixed 
on  all  edges. 


3.  Anisotropic  simply  supported  square  plate  under  uniform  normal  load. 


Conclusions 

The  boundary  residuals,  in  each  of  the  problems  studied,  exhibit 
a  steadily  decreasing  pattern  with  increasing  levels  of  truncation  of  the 
boundary  identities.  The  values  of  the  residuals  as  percentages  of  the 
maximum  values  of  the  corresponding  quantity  occurring  over  the  plate 
are  well  within  the  limits  to  which  the  boundary  conditions  can  practically 
be  specified.  The  very  rapid  convergence  of  the  fields  in  the  case  of  the 
cantilever  plate  may  be  attributed  to  the  use  of  a  symmetrised  representation, 
this  being  a  consistent  feature  of  the  use  of  the  symmetrisation  principle. 

In  the  case  of  the  fully  fixed  and  simply  supported  plates,  where  the 
symmetrised  representation  was  not  used,  the  symmetry  in  the  computed  fields 
becomes  more  exact  with  increasing  truncation  level.  The  consistent 
preference  shown  by  the  solver  routine  for  vertex  functions  as  opposed  to 
polar  functions  reinforces  the  expectations  of  accuracy  and  efficiency  for 
the  vertex  functions. 

The  relatively  small  number  of  degrees  of  freedom  required  to  obtain 
rapidly  convergent  solutions  demonstrates  the  efficiency  of  the  asymptotic 
field  representation  approach,  especially  in  the  case  of  the  severe  singular 
field  behaviour  exhibited  in  the  case  of  the  cantilever  plate  problem. 
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2. 


"The  Edge-Function  Method  for  Free  Vibrations  of  Isotropic  Plates" 
Sheehy  Richard.  Ph.D.  thesis  University  College  Cork  July  1979 


Abstract 


Previous  research  on  the  isotropic  plate  vibration  problem  is 
reviewed  briefly.  The  Edge-Function  Method  for  steady-state  vibrations  of 
polygonal  plates  is  developed*  and  asymtotic  expansions  of  the  solution 
in  the  neighbourhood  of  a  corner  are  obtained  under  general  boundary 
conditions  and  named  yertex  Functions. 

The  boundary  conditions  lead  to  a  system  of  homogenous  linear 
equations,  the  non-trivial  solution  of  which  yields  the  frequency  spectrum. 
A  brief  outline  of  the  computer  program,  VI BRAT,  as  developed  for  plate 
vibration  problems  is  given  and  is  included  as  an  appendix.  A  represent¬ 
ative  selection  of  results  obtained  for  a  variety  of  plate  geometries  is 
presented  and  a  critical  comparison  with  published  results  is  given. 
Hitherto  unsolved  problems  for  irregular  quadrilaterals  and  polygons 
are  presented  as  the  touchstone  with  which  to  evaluate  t.oe  Edge  Function 
Method  for  Free  Vibrations. 

Illustrative  Examples 

Exampl es  of 

(3)  Fully  Fixed  Rhombic  Plate 

(4)  Quadrilateral  Plate  with  Mixed  Boundary  Conditions 

(5)  Pentagonal  Plate  with  Mixed  Boundary  Conditions, 

together  with  Conclusions  are  attached  fron  the  above  thesis  to 
demonstrate  the  power  of  E.F.M.  in  Free  Vibration  problems. 

-  Fully  Fixed  Fdio-ribie  Flats: 
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Quadrilateral  with  Mixed  Boundary  Conditions 


(-1,0.5) 


(-1,-0. 5) 


(2,0.5) 


(1,-0. 5) 


FIG. 7. 5 


The  region  in  Fig. 7. 4  with  one  side  clamped  and  the  crr.er 
three  simply  supported  was  next  investigated.  Values  for 
frequencies  obtained  by  VIBRAT  and  presented  in  Table  7.21  are 
not  compared  with  other  investigators , since  a  search  of  the 
literature  did  not  unearth  such  results . 


Mode  and  )\  values 
~  j  ~  3~ 


E.F.M. 
i  L=3 


56.197  j  97.35  153.62; 


TA.BLB  7.21 


Percentage  R-M.S.  Boundary  Residuals  for  the  1  *"  Frequenc 
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MT 
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j 
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1 

1 

0.9768  x  lO-2 

1 

2 

0.9293  x  lO-2 

1 

3 

0.7200  x  lO-2 

1 

4 

0.8577  x  10“2 

2 

1 

0.2323  x  10-2 

3 

2 

0.1735  x  lO-1 

3 

3 

0.5584  x  lO  ^ 

3 

4 
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TABLE  7.24a. 
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.37 45 
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.1528 

3 
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4 
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1 

-5175 

2 
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5 
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Residuals 
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0.7341 
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8. 


CONCLUSION 


The  suite  of  programs  included  in  Appendix  B  are  written 
in  general  terms  and  allow  for  a  “solution  mix"  of  either 
scheme  A  -  Edge-functions  and  Fractional-Edge  functions  — 
or  scheme  B  -  Edge-functions  and  Vertex  functions. 

■ — 


5a- 

}  3  tec 


Each  solution  obtained  is  an  exact  solution  of  a 
ponding  problem  with  corresponding  boundary  values  as 
from  the  solution  sector,  a  measure  of  whose  deviation  irto. 
those  prescribed  is  provided  by  the  corresponding  root-rate 
square  boundary  residuals.  Accordingly, each  solution  provides 
a  Mathematical  Model  for  the  given  physical  problem.  The 
R-M.S.  values  provide  a  practical  criterion  for  the  engineer 
to  judge  the  acceptability  of  the  mathematical  model.  If  the 
mathematical  model  l:es  within  the  scatter  within  which  an 
engineer  can  specify  the  physical  problem,  then  the  mat  her. 
model  provides  as  “exact”  a  solution  to  the  physical  pros 
can  be  obtained-  A  scatter  of  tl%  is  used  in  this  thesis 


oical 

ers  as 


The  results  in  this  thesis  show: 

(1)  The  frequencies  with  either  "solution  mix'*  are  in  excellent, 
agreement  with  published  results.  In  cases  where  results  were 
not  available  there  is  excellent  agreement  between  values 
obtained  for  various  levels  of  harmonic  truncation. 

(2)  Boundary  Residuals  Rectangular  Regions:  ....... 

- - - - -  The  R.K.S. boundary 

residuals  for  rectangular  regions  were  quite  acceptable  using 

either  "solution  mix"  with  even  as  low  a  harmonic  truncation 

level  as  L-2. 

(3)  k -M^S.Residuals^^noi^RectanguIai^Reqionsj:  Scheme  A 

The  combination  cl  Edge  Functions  and  Fractional  Edge 
Functions  did  nor  produce  boundary  residuals  in  the  case  or  non- 
rectanguJar  regions  within  acceptable  limits  and, furthermore, 
the  required  symmetry  in  the  associated  mode  shapes  was  not 
satisfactory  1 y  reproduced . 
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(4)  BoundarY_RGSiduals_-_non_Rectan2ular Regions :  Scheme _B 

The  introduction  of  Vertex  functions  of  the  form  (4-13) 

.  reduced  the  boundary  residuals  by  a  factor  between  10  and  100 
from  those  obtained  using  Scheme  A,  and  accordingly  the  number 
t  of  harmonic  sets,  or  the  size  of  the  matrix,  required  to 

obtain  residuals  within  a  scatter  of  1  l%,can  be  correspondingly 
reduced. 

(5)  It  can  be  concluded, from  the  array  of  examples  presented,' 
that  Scheme  B,  the  Edge-Function  method  involving  r:;s  Vertex 
Function  developed  in  this  thesis,  is  both  satisfactory  and 
practical,  involving  relatively  moderate  computer  requirements. 
It  determines  acceptable  numerical  values  for  the  frequencies 
of  free  vibration  and  the  associated  mode  shapes  for  general 
polygonal  shaped  plates  with  mixed  boundary  conditions. 

f 

The  functions  in  Chapters  III  and  IV  may  appear  complicated 
j  but, in  keeping  with  Edge-Function  philosophy, the  computer 

formulae  (3-30)  and  (4-15)  provide  an  easy  means  of  controlling 
the  resulting  algebra.  Accordingly,  edge  functions  can  almost 
be  as  readily  applied  in  a  computer  program  as, say,  the 
trigonometric  functions.  The  touchstone  of  any  numerical  method 
must  be  the  ease  with  which  the  relevant  theory  can  be  pro¬ 
grammed  and  the  provision  of  a  practical  test  to  enable  an 
engineer  to  decide  whether  or  not  to  accept  the  results.  The 
Edge-Function  method  for  Vibration  problems  is  presented  for 
evaluation  on  the  above  criteria  and  especially  on  its  success 
in  solving  quadriloteral  and  polygonal  plate  problems,  as  in 
examples  (7-.'ll)  and  (7-2?.)  which,  to  the  author’s  knowledge, 
have  not  been  solved  hitherto. 
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"Analysis  of  Elastic  Plates  using  the  Edge-Function  Method" 
Ahmed  A. A.,  Ph.D.  thesis  Univ.  of  Dundee  October  1979 

Abstract 


This  thesis  describes  an  investigation  into  the  applicability  of 
the  edge-function  method  for  the  engineering  analysis  of  plates  and  slab 
bridge  decks.  The  specific  objective  is  the  assessment  of  accuracy  .. 
attainable  economically  in  evaluating  deflections,  bending  moments  and 
shears  for  the  complex  situations  encountered  by  the  bridge  designers. 

First,  the  development  of  the  edge-function  theory  is  presented 
for  the  analysis  of  thin  plates  and  slabs  with  different  planforms, 
types  of  loading  and  support  conditions.  The  bending  problem  is  formu¬ 
lated  for  both  isotropic  and  orthotropic  plates,  and  general  solution 
forms  are  obtained.  The  technique  is  associated  with  solutions  from  all 
critical  regions  of  the  boundary  and  the  presence  of  corner  functions  is 
found  necessary  to  cater  for  singularities  and  to  accelerate  the  conver¬ 
gence. 


Two  computer  programs  were  developed  to  facilitate  analysis  of 
isotropic  and  orthotropic  plates  and  slab  bridge  decks;  they  provide  for 
the  computation  of  column  reactions  as  well  as  bending  and  twisting 
moments  at  any  required  station.  A  brief  outline  of  the  computer  programs 
is  given. 

Numerical  examples  are  given  along  with  extensive  comparisons 
with  the  results  of  other  authors  using  other  numerical  methods  and  tests 
on  models.  The  investigation  shows  excellent  agreement  and  indicates 
that  a  generally  accurate  and  inexpensive  method  of  solution  is  now 
available.  The  advantages  of  this  technique  over  the  conventional  domain 
type  methods  (finite  element  method  and  finite  difference  method),  which 
require  division  of  the  whole  plate  into  elements,  lie  in  the  much 
smaller  size  of  the  matrix  to  be  solved,  and  it  is  thus  more  economical 
and  suited  for  programming  on  medium  size  computers. 

Recommendations  for  further  development  and  application  of  the 
method  to  ther  types  of  bridge  deck  are  stated. 


“The  Edge-Function  Method  for  Elliptical  Cracks  in  a  Prismoidai  Body* ‘ 
Quinlan  P.M.,  Grannell  J.J.,  Atlufi  S.N.,  Fitzgerald  J.E. 

Proceedings  Second  International  Symposium  on  Innovative  Numerical 
Methods  in  Engineering  Science  Montreal  (1980). 

Abstract 

Stress  fields  are  developed  to  model  the  behaviour  of  Elliptical 
cracks  in  a  prismoidai  body  using  confocal  potential  functions. 

Stresses  and  displacements  are  obtained  in  a  form  suitable  for  comput¬ 
ation  without  having  to  introduce  Jacobian  Elliptic  Function. 

A  special  quadrature  subroutine  has  been  developed  for  Crack 
Functions  and  they  are  now  incorporated  into  computer  program  PQDISK  [5] 
for  "3-Dimensional  Stress  Analysis  using  The  Edge-Function  Method". 

Applications  are  made  to  a  wide  range  of  buried  elliptic  cracks, 
both  large  and  small,  under  concentrated  and  distributed  normal  loadings* 
The  effects  of  increasing  the  number  of  Crack  Functions,  and  the  corresp¬ 
onding  boundary  residuals,  are  studies  and  several  illustrative  example 
are  given. 


INTRODUCTION 

Since  the  publication  by  Shah  and  Kohayashi  [1]  in  1971  which 
developed  the  confocal  potential  functions jintroduced  by  Segedin  [2] 
for  potential  problems,  to  determine  the  Stress  Intensity  Factor  for 
an  Elliptical  Crack'  in  an  infinite  body  under  arbitrary  normal  loading 
the  extension  of  that  work  to. finite  bodies  has  presented  a  major 
challenge.  A  boundary  integral  solution  was  given  by  Cruse  [3]  in 
1975  which  v/as  followed  in  1977  by  the  hybrid  element  solution  of 
Kathiresan  [4]  with  a  combination  of  finite  and  special  crack  tip 
el emonts . 

However  little  progress  appears  to  have  been  made  in  analytical 

methods  using  [1],  This  is  due  mainly  to  the  almost  prohibitive 

algebraic  difficulties  in  expressing  in  terms  of  Jacobian  elliptic 

functions  all  the  second  and  third  derivatives  of  the  stress  potentials 

<j>..  reouired  for  the  resulting  boundary  stresses »' and  difficulties  in 

if 

evaluation.  Paper  [13  obtains  but  one  derivative,  6  ^ij  ,  and  that 
for  only  i  +  j  s  3. 

The  present  paper  presents  a  straightforward  calculation  scheme 
for  all  the  required  derivatives  of  6..  without  introducing  any 

•  u 

Jacobian  elliptic  functions  and  no  particular  difficulties  are 
encountered  with  values  of  i  and  j  in  the  range  (0,  10).  The  resulting 
Elliptic  Crack  Functions  have  beer,  included  in  the  author's  computer 
program  PODISK  15]  for  Three  Dimensional  Stress  Analysis  in  Prismoidal 
Bodies  and  numerical  results  are  given  for  both  small  and  large  cracks. 
In  sharp  contrast  to  [33  and  [4]  large  cracks  can  be  handled  as  readily 
as  small  cracks.  The  computer  times  required  are  at  least  an  order  of 
mangitude  less  than  those  in  [33  and  [43;  stresses  and  displacements 
are  readily  calculated  at  any  point,  and  each  additional  load  case 
reouires  but  a  few  percent  of  the  time  for  the  first  case  when  all- 
cases  are  processed  together. 


A  Crack  function  D .3] capability  has  been  inserted  in  computer 
program  PODISK  £5]  which  implements  "The  Edge  Function  Method  for 
Prismoidal  Bodies"  and  any  required  number  from  1  to  15  of  such 
functions,  can  be  included  in  the  solution  scheme  for  modelling  an 
elliptical  crack. 

The  Edge  Function  Method  is  described  in  16,  7,  8,  9,  11,  12],  and 
space  prevents  any  further  elaboration  here  except  to  point  out  that 
each  boundary  condition  is  imposed  on  each  face  in  either  a  continuous 
or  a  discrete  least  squares  sense.  Consequently  the  boundary  condition 
of  zero  normal  stress  in  the  cracked  ellipse  is  imposed  in  a  least 
squares  sense;  the  squares  of  the  residuals  in  the  cracked  ellipse 
being  minimised  w  r  t  the  coefficients  of  the  crack  functions.  The 
zero  shear  condition  is  provided  by  the  built  in  symmetries  of  the 
functions. 

A  cube,  2x2x2,  with  a  centrally  located  elliptical  plane  crack 
loaded  normally  on  the  faces  parallel  to  the  crack  plane  as  shown  in 
Figure  2,  is  considered.  Four  different  load  types  -  (!)  concentrated 
load  (2)  line  load  (3)  natch  load  modelled  by  equivalent  concentrated 


loads  using  the  Bousinesque  solution  and  (4)  a  linear  stress  field- 
producing  a  uniform  load^can  be  taken  together  in  each  job;  the 
four  loading  cases  processed  together  taking  only  about  15%  more 
computer  time  than  would  be  reouired  to  obtain  a  single  solution  on 
its  own. 

Four  different  problem  sets  are  presented  under  to  illustrate 
the  use  of  the  crack  functions  - 


(a)  Penny  Shaped  Crack 

As  given  in  Green  and  Sneddon  UO]  the  fracture 
coefficient  Kf  for  a  penny  shaped  crack  of  radius 
a  in  an  infinite  solid  in  a  uniform  stress  field  p 
normal  to  the  plane  of  the  crack  is  -  0 

Kj  =  2p0/a77  (7.1) 

The  fracture  coefficient-  K.  was  obtained  for  several 
values  of  the  radius  a  and  the  ratio  K,/K^,  which 
measures  the  boundary  effects  on  the  fracture  coefficient 
is*  given  under*  together  with  the  root  mean  square 
of  the  residual  stress,  Ap,  remaining  on  the  crack 
ellipse 


a 

VKi 

Ap 

.2 

1,02 

0,1 4*1 0"5 

.4 

1.11 

0.68*1 O'4 

.6 

1.39 

0.1 3*1 0~2 

.8 

2.34 

0.40*10"2 

• 

Table  2 

(b)  Small  Elliptical  Crack 

The  problems  solved  by  Cruse  [3^Jand  Katiresan  M 
obtained  K,  for  relatively  small  ellipses  where 
ab  =  .04  in  Figure  2  with  a,;  =  1,  2,  3,  4  and 
under  a  uniform  strees  fielB.  Program  PODISK  was 
run  for  these  cases.  -Pesults  are  given  in  Table  3 
under  for  Kj/K|  at  10  points  on  a  quadrant  of  the 


;  1 

*  l 


vp  --'<1  -  ■ 
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crack  front  corresponding  to  equal  intervals  in 
the  elliptical  angle  o  for  o(0,7iy2),  an(*  f°r 

residual  stress  Ap  on  the  crack  ellipse.  As  given 
in  Kathiresan  [4],  analogous  to  (7.1),  K?  for  all 
ellipse  of  semi-axes  a  and  b  is  given  by- 

K*  =  P0*^  (sin2e  +  ,b.2  cos2e>^  ,  a  >  b  p 

w  a 


# 

E(K)  J 


'1  -  k  sin  k  =1 


2=1  -  bU 


Ratio  K-j/K,  at  10  points  on  Crack  Front 


1.02  Constant  at  all  points  0.14  *  10~^ 
.932  .952  1.00  1.07  1.14  1.20  1.25  1.29  1.31  1.32  0.13  *  10“5 
.828  .  874  .  976  1.09  1.19  1.28  1.35  1.39  1.42  1.43  0.23  *  10~5 
.746  .819  .961  1.11  1.22  1.32  1.39  1.45  1.48  1.49  0.14  *  10~4 


Table  3 


The  results  all  agree  to  within  2%  with  those  given 
by  Kathiresan  [4],  They  were  obtained  using  6  crack 
functions,  together  with  60  plane  functions  to  model 
the  boundary  effects  and  a  linear  stress  field.  The 
computer  time  reouired  for  e3ch  ellipse  is  approx¬ 
imately  six  times  that  reouirad  to  solve  70 
simultaneous  ecu at ions  by  ordinary  Gaussian 
elimination  and  consequently  is  an  order  of  magnitude 
less  than  was  reouired  in  13]  or  [4],  The  eouations 
took  about  five  times  as  long  to  generate  as  to  solve. 

(c)  Elongated  Elliptical  Crack 

As  a  demonstration  of  the  power  of  program  PODISK 
examples  (c)  and  (d)  are  presented.  The  first 
takes  a  very  elongated  crack  with  a  =  0.8  and 
b  =  0.05.  It  uses  120  plane  functions  to  model 
the  boundaries  ar.d  10  crack  functions  derived 
from  D. 3] with  the  even  values  for  i  and  j  as 
under  to  provide  the  necessary  symmetries  - 
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(7.3) 
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Solutions  are  readily  obtained,  in  the  sane  job,  using 
N  ~  1,  2,  3  ...  10  crack  functions  and  for  the  four 
load  types  mentioned  above.  Table  4  shows  the  variations 
of  the  residual  stress  Ap  cn  the  crack  plane  for  N  crack 
functions  and  each  load  type,  while  Table  5  gives  some 
illustrative  K,  factors  at  points  1,  3,  5,.  7,  9  and  10 
as  in  Figure  3,  on  a  quadrant  of  the  ellipse  for 
concentrated  and  patch  loads. 


Concen.  Load  at  (.2,  .3) 

Line  Load  from 
(0,.4)  to  (.5,. 4) 

Patch  Load 
on  1x1 

Uniform  Stress 

.4110  00 

.1720  00 

.151D-01 

.2750-03 

.4070  00 

.168D  00 

.1500-01 

.2720-03 

.2690  00 

.1710-01 

.505D-03 

.1840-4)4 

.2660  00 

.169D-01 

.505D-03 

.1820-04 

.2660  00 

;169D-01 

.5060-03 

.1810-04 

.1120-01 

.438D-02 

.3620-03 

.1250-04 

.1120-01 

.4360-02 

.362D-03 

.1250-04 

.1110-01 

.434D-02 

[ 

.3580-03 

.1200-04 

.110D-01 

.434D-02 

.3520-03 

.1170-04 

.9210-02 

.4330-02 

.2710-03 

.7490-05 

Table  4  Stress  Residuals  on  a  small  Crack  Plane 
(Effect  of  increasing  the. number  of 
Crack  Functions! 


Concentrated  Load 

K-j  Factors  at  Points  1,  3, 

5,  7, 

9,  10 

Nc  1 

.843 

1.99 

2.71 

3.14 

3.35 

3.37 

2 

.842  • 

1.98 

2.68 

3.08 

3.26 

3.29 

3 

.351 

.994 

1.92 

2.98 

3.71 

3.81 

6 

.392 

1.04 

1.89 

2.93 

3.72 

3.84 

9 

.392 

1.04 

1.89 

2.93 

3.72 

3.83 

Patch  Load 
(Bousinesque) 

Nc  1 

.399 

.939 

1.28 

1.48 

1.58 

1.59 

2 

.398 

.939 

1.28 

1.48 

1.58 

1.59 

3 

.380 

.902 

1.25 

1.48 

1.60 

1.61 

5 

.380 

.902 

1.25 

1.48 

1.60 

1.61 

Table  5  Ki  Fracture  Coefficients 

(Crack  Front  Oistritution) 
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Table  'l  shows  that  six  crack  functions  are  required 
to  reduce  the  residual  stresses  to  Iff-  in  the  . 
concentrated  load  case,  while  three  functions  reduce 
the  residuals  to  less  than  .01?=  in  the  case  of  a 
patch  load.  The  variations  in  the  K,  factors 
correspond  to  the  residual  stresses  and  no 
significant  change  occurs  in  the  values  after  six 
and  three  crack  functions  respectively.  The 
residual  stresses  on  the  surface  of  the  cube  were 
less  than  of  the  applied  loading  and  could  be 
reduced  further  by  increasing  the  number  of  plane 
functions  used  in  the  modelling. 

(d)  Large  Elliptical  Crack 

Program  P00ISK  v/as  run  for  a  crack  with  a  =  0.8, 
b  =  0.7  usino  the  same  number  of  functions  as  in 
(c)  and  results  analogous  to  those  in  Table  4 
are  given  in  Table  6. 


Cone.  Load 
at  (.2,  .3) 

Line  Load  from 
(0,  .4)  to  (.4,  .4) 

Patch  Load 
on  1  x  1 

.6560  00 

.2230  00 

.8920-01 

.3170  00 

.1690  00 

.454D-01 

.4830-01 

. 1700-01 

.272D-02 

.248D-01 

.1600-01 

.224D-02 

.1210-01 

.135D-01 

.2030-02 

.6130-02 

.2570-02 

.120D-02 

.5500-02 

.1890-02 

.113D-02 

.543D-02 

.1890-02 

.1130-02 

.5490-02 

.170D-02 

.113D-02 

.4740-02 

.134D-02 

.5890-03 

Table  6  Stress  Residuals  on  a  Large  Crack  Plane 
(Effect  of  increasing  the  number  of 
Crack  Functions) 


The  K-j  factors  corresponded  to  the  residual  stress 
•pattern  in  a  similar  manner  to  the  factors  in  Table 
5.  It  is  evident  that  five  and  three  crack  functions 
are  required  for  concentrated  and  patch  loads 
respectively  to  give  K-j  factors  to  within  It. 

Work  is  progressing  rapidly  in  adapting  PQDISK  to 
deal  with  surface  and  corners  flaws  and  also  to 
deal  with  fracture  of  lfode  II  and  III  since  the 
corresponding  crack  functions  can  be  generated 
from  confocal  harmonics  (1.3).  Results  are 
expected  in  time  for  the  Montreal  meeting. 

Displacements  and  stresses  are  readily  calculated  at 
any  point  in  the  body,  and  they  are  consistent  with 
the  boundary  residuals  reported  in  this  section. 

Space  limitations  prevents  their  presentation  in  this 
paper. 
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